Proceedings: Evolution of Large Scale Structure - Garching, August 1998 



CONSTRAINTS ON COSMOLOGICAL PARAMETERS 
FROM CURRENT CMB DATA 



James G. Bartlett, Alain Blanchard, Marian Douspis, Morvan Le Dour 
Observatoire de Strasbourg, 61000 Strasbourg, FRANCE 



ABSTRACT. We discuss the constraints one can place on cosmological pa- 
rameters using current cosmic microwave background data. A standard x 2 - 
minimization over band-power estimates is first presented, followed by a dis- 
cussion of the more correct likelihood approach. We propose an approxima- 
tion to the complete likelihood function of an arbitrary experiment requiring 
only limited and easily found information about the observations. Examina- 
tion of both open models - (H, h,Q,n) — and flat models (Q + SI a = 1) ~~ 
(O, r2j,, h, Q, n) — leaves one rather robust result: models with small curvature 
are favored. 



1 Introduction 

Current cosmic microwave back- 
ground (CMB) results already per- 
mit some constraints to be placed 
on certain cosmological parameters, 
such as the density parameter, £1, 
the Hubble constant, H = ?i(100 
km/s/Mpc), etc... (Bond & Jaffe 
1996; Lineweaver et al. 1997; Han- 
cock et al. 1998; Bartlett et al. 
1998a). The present data set is most 
clearly summarized in the power 
spectrum plane as a set of (flat) 
band-power estimates, as shown in 
Figure 1. Presence of the "Dopplcr 
peak" is indicated by several obser- 
vations on its upward slope, includ- 
ing the new QMAP results (Oliveira- 
Costa et al. 1998), and most impor- 
tantly the Saskatoon points (Netter- 
field et al. 1997). It is worth em- 
phasizing, all the same, that MSAM 
(Cheng et al. 1997) is not as support- 
ive of the peak as Saskatoon, and 
these are the only two experiments 
covering the peak at the present 
time. Analysis of newly obtained 
data, e.g., from BOOMERANG'S 
North American flight, MAXIMA or 
MAT, all capable of measuring the 
peak with good precision, will help 
to resolve the issue. 

We have used these power spec- 
trum estimates to constrain cos- 
mological parameters in both open 
(fl A = 0) and flat (C + fl A = 1) 
scenarios within the context of infla- 
tionary models with cold dark mat- 
ter. After presenting some of our re- 
sults based on a x 2 ~ mm i m i za ti° n j 
we discuss the shortcomings of this 



approach. An easy-to-use approxi- 
mation to the full likelihood func- 
tion of an experiment is then pro- 
posed and appears to lead to less re- 
strictive constraints than indicated 
by the \ 2 — minimization. Neverthe- 
less, one result remains the same, de- 
spite the general differences between 
the two methods: purely open models 
with Q < 0.4 are strongly disfavored, 
arguing against a Universe with large 
curvature. 

2 x 2 ~ approach 

The simplest 
approach is a straight-forward appli- 
cation of the x 2 ~statistic to the data 
shown in Figure 1. We used CMB- 
FAST (Zaldarriaga et al. 1998; Sel- 
jak & Zaldarriaga 1996) to construct 
a grid of power spectra spanned by 
(Q,h, Q,n), for purely open mod- 
els, and by (Q, {If,, h, Q, n) for flat 
models with Q + Q\ = 1. In Fig- 
ures 2 and 3 we show the resulting 
constraints on Q and h as dashed 
lines for both types of models. The 
confidence regions are defined rela- 
tive to the minimum value of \ 2 ; 
for example, over two parameters the 
68% and 95% contours are defined by 
the projection onto the plane of the 
multidimensional ellipse inscribed by 
xLn + 2 -3and X 2 min + 6-17 respec- 
tively. at The most striking result is 
the lower limit on the density param- 
eter in open models - Q > 0.5 at 
95% confidence. This is quite clearly 
due to the apparent position of the 
Doppler peak around I = 300 in the 
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Figure 1. Current CMB power spectrum estimates. The curve shows the 
predictions of a model with fl = 0.9, h = 0.35, Ui,h 2 = 0.015, Q = 17^K and 
n = 0.94. 
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Figure 2. Open model constraints. The dashed lines show the x 2 - 
minimization constraints (68% and 95% confidence boundarys), while the 
solid lines correspond to the approximate likelihood function. 



data. In this sense, the current data 
do not favor large spatial curvature, 
which would move the peak too far 
to the right in the figure. An illustra- 
tive example of this is given in Fig- 
ure 4, where we see how miserably a 
somewhat "conservative" model with 
Q = 0.2 and h = 0.6 performs. 

3 Of Likelihood Functions 

The common method just em- 
ployed does not strictly apply to the 
CMB data in Figure 1, for it sup- 
poses that the data points are gaus- 



sian distributed, and this is certainly 
not the case. Even if the tempera- 
ture of individual pixels on the sky 
can be considered a gaussian ran- 
dom variable, as predicted by infla- 
tionary scenarios, if the experimen- 
tal noise is also gaussian, a band- 
power estimate will not be gaussian, 
because it represents the variance of 
the pixel fluctuations. The variance 
of a gaussian random variable is not 
itself gaussian distributed. Thus, at a 
fundamental level, the x 2 ~ a PP roac h 
is not suited to our problem. 

We must employ a more general 
likelihood function to properly con- 
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Flat models 




Figure 3. Flat model constraints. The lines have the same meaning as in 
Figure 2. 



strain our parameters. This may be 
done for a given experiment if one as- 
sumes that the temperature fluctua- 
tions on the sky, the pixel values, are 
indeed gaussian random variables. In 
this case, with a set of N plx pixel val- 
ues arranged in a data vector d, we 
may write the likelihood function for 
the parameters, represented by the 
vector = {Q, h, ...}, as 



£(G) = Prob(d\B) 



( 27r )JVpix/2| C |l/2 



The key object is the covariance ma- 
trix of the pixels, C: 



CV, = < didj >= C(6» i3 ) + Nij (2) 
i 

xPticosOij) + 



where Wi is the window function de- 
fined by the experimental beam and 
Nij is the experimental noise covari- 
ance matrix. We have imagined a sit- 
uation were the observations are rep- 
resented by a set of simple pixel val- 
ues, e.g., a map, but the argument 
is the same for temperature differ- 
ences, or any linear combination of 
pixel values. 

The correct manner to proceed 
would be to obtain C for each exper- 
iment (including all correlations be- 
tween experiments covering the same 



sky zone) and form the complete like- 
lihood function for the entire data 
set. This is rather demanding and re- 
quires sometimes hard to find exper- 
imental information, such as a full 
noise covariance matrix. We have de- 
veloped a general purpose approx- 
imation to the likelihood function 
which, by comparison to the full like- 
lihood functions from COBE, Saska- 
toon and MAX, seems to work well 
j-^-j (work in progress, and reported in 
Bartlett et al. 1998b; for other work, 
see Bond et al. 1998 and Wandclt 
et al. 1998). The approximation re- 
quires only limited information con- 
cerning a given experiment, such as 
the number of pixels, the band- 
power estimate and error bars, and 
an idea of the noise level — all readily 
found in the literature. 



4 Results from the Approxi- 
mate Likelihood 

Applying our approximate likeli- 
hood ansatz to each data point in 
Figure 1, we have revisited the pa- 
rameter constraints imposed by the 
data set. The results in the (H, h)- 
planc arc shown in Figures 2 and 3 as 
the solid lines; once again, these are 
projections of ellipses defined relative 
to the maximum value of the likeli- 
hood. We see that in general the con- 
fidence regions are larger than those 
deduced from the x 1 — approach. This 
is most spectacular for the flat mod- 
els, where virtually all of the parame- 
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Figure 4. A model with Q = 0.2, H a = 60 km/s/Mpc, n = 1, Q = 20 fj,K, 
fl b h 2 = 0.015 



ter space in this plane is now allowed 
at "95%" confidence. One should be 
cautious interpreting these results, 
however, for a more correct definition 
of confidence regions would integrate 
the likelihood function over the pa- 
rameters not shown. In any case, the 
difference in the contours, defined in 
the same way, clearly indicates that 
the approximate likelihood function 
has a different distribution over the 
parameter space than that suggested 
by the x 2 — approach. 

Thus, we see that care must be 
taken in drawing conclusions from a 
X 2 — minimization to band-power es- 
timates. Nevertheless, the best fit 
models do not change drastically. We 
also find about the same lower limit 
as before to Q in open models - Q > 
0.4. 



5 Conclusion 

It is important to realize that the 
CMB is already beginning to di- 
vulge its treasures of information. 
The method used to extract quan- 
titative constraints on cosmological 
parameters requires some care, and 
a simple \ 2 ~ analysis using band- 
power estimates is not the best ap- 
proach, leading to constraints which 
appear tighter than the data war- 
rant. One important result already 
indicated by the present data is that 
the spatial curvature of the Universe 
cannot be too large; this is quantified 
by the lower limit of Q > 0.4 for the 



open models examined here. 
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